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Abstract 

We consider 3-point numerical schemes for scalar Conservation Laws, that are oscil- 
latory either to their dispersive or anti-diffusive nature. Oscillations are responsible for 
the increase of the Total Variation (TV); a bound on which is crucial for the stability of 
the numerical scheme. It has been noticed ( [AKMOlj . [AMT04j . [AMSOSj l that the use 
of non-uniform adaptively redefined meshes, that take into account the geometry of the 
numerical solution itself, is capable of taming oscillations; hence improving the stability 
properties of the numerical schemes. 

In this work we provide a model for studying the evolution of the extremes over non- 
uniform adaptively redefined meshes. Based on this model we prove that proper mesh 
reconstruction is able to control the oscillations; we provide bounds for the Total Variation 
(TV) of the numerical solution. We moreover prove under more strict assumptions that the 
increase of the TV -due to the oscillatory behaviour of the numerical schemes- decreases 
with time; hence proving that the overall scheme is TV Increase-Decreasing (TVTD). 

1 Outline 

We study the scalar Conservation Law in one space variable 

ut + f{u), = 0, xe [0,1], te[0,T] 

where the flux function / is considered to be smooth and convex. The initial condition through- 
out this work is: 

uo{x) = < ^' ^ - ^° , for G (0, 1) 
10, X > Xq 

We discretise spatially over a non-uniform, adaptively redefined mesh, 

M'' ^ {0 = x'^ <x'i < ■■■ Kx"}^ = 1} 

The manipulation of the non-uniform mesh and the time evolution are combined into the 
Main Adaptive Scheme (MAS): 

Definition (Main Adaptive Scheme (MAS)). Given, at time step n, the mesh A/" = {a = 
x" < • • • < x^^ = b} and the approximations [/" ~ {^^i j ■ • ■ j ^^n}' steps of the (MAS) are as 
follows: 
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1. (Mesh Reconstruction) 

Construct new mesh Af"+^ = {a = < • • • < x^"^ = h] 

2. (Solution Update) 

Using the old mesh M" the approximations t/" and the new mesh M""*"^; 

2a. construct a piecewise linear function V"'{x) such that V^"|j\fn = U" 
2b. define the updated approximations U" — {w", . . . jU^} as U" = 

3. (Time Evolution) 

Use the new mesh M"+^, the new approximations [/" and the numerical scheme to march 
in time and compute U^^^ = {u"^^^ , . . . , u^^} 

4. (Loop) 

Repeat the Step 1.-3. with M^~^-^, [/"+^ as initial data. 

We shall discuss the MAS in more details in section 21 for the moment we shall make some 
brief notes on the Steps 1.-3. 

Remark 1. Regarding the mesh reconstruction (Step 1.), we note that it is performed in each 
time step and is the crucial ingredient of the MAS. This procedure creates a new non-uniform 
mesh A'/"+"' with the same number of nodes as the previous one M" . The construction of the new 
mesh depends on the geometry of the numerical solution itself. Regarding the solution update 
(Step 2.) and time evolution (Step 3.), we note that in this work we consider interpolation over 
piecewise linear functions and Finite Difference schemes. The proper setting for conservative 
solution update and Finite Volume schemes will be addressed in a different work. 

The objective of this work is to place conditions on the steps of MAS such that the resulting 
numerical solutions are TV stable even when oscillatory numerical schemes are used for the 
time evolution (Step 3.). More specifically, we shall prove that proper non-uniform meshes are 
able to tame the TV increase due to oscillations, furthermore we shall prove that in some cases 
the TV increase due to oscillation decreases with time; hence yielding a Total Variation Increase 
diminishing (TVID) scheme. 

In section [5] of this work we list and explain the requirements that we place on the MAS. In 
section [31 we discuss the creation and propagation of oscillations at the level of extremes. We 
present a model for the extremes that takes into account the steps of MAS. Based then on the 
model we prove the TV results of this work. In section [H we discuss the Main Adaptive Scheme 
(MAS) in more detail. We analyse the way non-uniform meshes are constructed and how the 
numerical solution is updated over the new mesh. We discuss the numerical implementation of 
the MAS and present some graphs depicting its basic properties. In section [5] we discuss the 
numerical implementation of the requirements introduced in section [51 In section [H] we perform 
numerical tests, where we consider known oscillatory numerical schemes and prove that these 
schemes satisfy the requirements introduced in section [51 We provide comparative numerical 
results for both uniform and non- uniform meshes. 

2 Requirements 

In this section we present the requirements that we place on the steps of the MAS. We once 
again mention that we perform interpolation over piecewise linear functions for the Solution 
Update (Step 2.) and that we use oscillatory Finite Difference schemes for the Time Evolution 
(Step 3.). The proper setting for dealing with Finite Volume schemes with a conservative 
reconstruction shares many things in common with this work but will be presented separately. 
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Let M" = {x^,i = 1, . . . , N} be the initial mesh, [/" = {w", . . . , u^} the initial approxi- 
mations at the time step n, and = {^j^^ij = 1, • • • , N}, [/" = {w", . . . , u^} be the new 
mesh and updated approximations yielding from the Steps 1. and 2. of the MAS. 

To introduce the first requirement, we recall that the development of this work docs not 
assume the use of any specific numerical scheme for time evolution (Step 3.). In the contrary 
the discussion that will take place and the proofs that will be presented arc valid for every 
numerical scheme that satisfies the Evolution requirement: 

Requirement 1 (Evolution requirement). There exists a constant C > independent of the 
time step n and the node i such that, 

|<+^-<| <Cmax{|<+i-^rU<-<-il} (1) 

Remark 2. The meaning of this requirement is that the numerical scheme, which is responsible 
for the time evolution (Step 3.), does not produce abrupt results. We shall see in the examples 
addressed in section that the constant C is an increasing function of the CFL condition. For 
every scheme that we will use, we shall prove this requirement and also compute the dependence 
of the constant C to the CFL condition. 

We move on to the second requirement, which is placed on the mesh reconstruction procedure 
(Step 1). We first start with a definition. 

Definition 2.1. We say that the approximate solution [/" = {wf, « = 1, . . . , A'^} defined over 
the mesh A/" = {2^";* = 1, . . . , A} exhibits a local extreme at the node cc" if > 
(local maximum) or < uf_i,uf_^i (local minimum). 

Requirement 2 (Mesh requirement: A-rule). There exists a constant < A < 1 independent 
of n and i such that, if x^'^^ £ [a;f , and U"" exhibits a local extreme at the node i (resp. 

i -\- \) then 

x]+' - x2 > (I - A)(.tIVi - (^respecttvely x^^, - x]+' > (1 - X){x^+, - x^)) (2) 

or equivalently 

xr+i - xj+i < A(a;r+i - (respectively x]+' ~ x^ < Xix^^, ~ x")) (3) 

Remark 3. The meaning of the X-rule requirement, Rel. is that the new nodes x^'^^ avoid 
the places of the old extremes x" ( or x^^_^-^) by at least 1 — A 0/ the respective interval (x", xf^i). 

The A-rule Requirement is placed on the mesh reconstruction but also affects the values 
of piecewise linear functions, we call this A-rule effect. The following remark discusses their 
relation. 

Remark 4 (A-rule effect for piecewise linears and interpolation). Assume that u(x) is a piece- 
wise linear function that oscillates as depicted in Fig.^. Assume moreover that the new nodes 
respect the X-rule Req. 0) at the extremes in the respective subintervals. Let also y = v be a 
horizontal line that separates the extremes. According to Fig.(^, x^^l G [2;f_i,a^"] and by the 
X-rule requirement 

x]^l-xU>il-X)ix^-xti) 

we get 







X X 


□ X X 


+ (1 - A)Aa;? 


- (1 - A)Axf 



Figure 1: The A-rule states that the new node x"~^^ £ should avoid the old node 

x" (respectively xf^^), where the approximate solution might exhibit an extreme, in the sense: 
xf + {1- X)Axf < a;"+^ (respectively < xf+i - (1 - A)Ax^), where Axf = xf^^ - xf 

Since u is linear in the interval [x'^_i,xf] 

uix-)~u{x]+l) _ x^-x]+l ^ 
u(a;f) -u(a;f_i) - a;f_i 

by the monotonicity of u in the interval [xf_i,xf] the previous relation recasts into 

u{x]+l)^u(x-)<\(u{x-_,)~u{x-)) 

since < A < 1 and it(a;") < v the previous relation reads 

uix]+l) < Au(a;,r_i) + (1 - A)u(a;^) < Am(x^_i) + (1 - A)w => u{x]+l) -v< A(u(x^„i) - v) 

For the rest of the extremes of Fig. that is u{x''l) and u[x^j^i) and for the respective new 
nodes x^^^ , x^^l we can similarly prove that, 

Hx']^^) - v\ < \\u{x'l) - v\ and \u{x';Xl) ~ v\ < \\u{x1+^) - v\ 

The meaning of this remark is the following: if a new node respects the X-rule in a specific 
interval, the same is true for the value of this new node with respect to the values of the function 
at the end points of this interval. 

The case where several nodes are placed between consequent extremes will be addressed in 
Section O There 

The last remark connects the A-rule -which is applied on the mesh- with the values of the 
piecewise linear function. In the rest of this work we shall refer to this connection as the X-rule 
effect. 

To pass to the third requirement, we first note that the overall phenomenon consists of two 
major steps, the mesh reconstruction (Step 1.) and the time evolution (Step 3.). We need 
to study these steps both separately and together. For the separate analysis the requirements 
Req.((T]) and Req.([2|) are sufficient, but for the joined analysis one more requirement is needed. 
The necessity for the extra requirement is explained in the following remark. 

Remark 5. The effect of time evolution (Step 3.) which is due to the numerical scheme, 
can be studied with several means, such as the respective Modified Equation of the scheme. In 
the contrary, the effect of the mesh reconstruction (Step 1.) cannot be analysed with classical 
methods, since it takes place between two consequent time steps. In other words, the mesh 
reconstruction is not related to the time evolution of the numerical solution. 
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Figure 2: This figure depicts the application of the A-rule in the case of a piecewise hnear 
function. The places of the new nodes are depicted along with the old extremes. 

The major contribution of this work is the merging of the effects that these steps have on the 
appearance and the evolution of oscillations; hence on the TV of the numerical approximation. 
The merging of these effects is quantified in the Coupling requirement: 

Requirement 3 (Coupling requirement). The constants C of the Evolution Req.^ and A of 
the X-rule Req.(^ are connected via the following relation: 

X + 3XC <1 (4) 

3 Time step analysis 

In this section we discuss the appearance and evolution of local extremes. We devise recursive, 
with respect to the time step n, relations regarding the evolution of the extremes based on the 
requirements presented in the previous section. We present and prove the theoretical results of 
this work; these include bounds on the extremes, bounds on the TV increase due to oscillations, 
we moreover prove that in some cases the TV increase, decreases with time steps n. 

3.1 Recursive relations 

The discussion regarding the creation and evolution of the extremes is performed in a time step 
by time step manner. In every time step we shall discuss their temporal evolution and their 
spatial modification. 

We start with a jump initial condition. Fig. ([3]), which wc discrctize over a non- uniform 
mesh. In the description that follows we have split every step into two sub-steps. The first 
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Figure 3: This is the initial condition. In this configuration we set ai = |?},° — 

sub-step is the time evolution, which is due to the numerical scheme and is governed by the 
evolution Rcq.(IT]) and the second sub-step is the spatial modification, which is due to the mesh 
relocation and the solution update procedure and is governed by the A-rule Rcq.Q and the 
couphng Req.([3]). 

1-st step We refer to Fig. ([31) for a graphical description of the configuration. The first nodal point 
located at the top of the shock, that is will evolve according to the Evolution Req.(IT]), 
which reads 

\ul-u^\<CmB.^{\u^-ul,l\u^~ul,\} 

Since we consider jump initial conditions, it is obvious that 

|h°-CiI=0 and \u^-ul,\<TV{u'>) 

We denote oi = |u° — "i+il and define ai ~ Cai, so the evolution Rcq.([T]) for the value 
M° reads, 

\ul - u',\ < Cmax{|zlO - Cil, I"? " < C\u^ ~ = Ca^ = a, 

1/2 

To introduce the notation for the continuation of this work we define to be the 

maximum magnitude of this extreme, hence 

El" = \u\-u'l\=a^ 

1/2 

To explain the symbolism, we use the letter E in E-^ because we refer to the magnitude 
of extremes, the subscript i states that we refer to 1-st extreme and the superscript ^/■^ 
states that we have moved from the time step fc = 0, with the use of the numerical scheme 
(time evolution step) but the mesh reconstruction step has not taken place yet. 

We refer to Fig. ([4]) for a graphical description. We then perform the mesh reconstruction 
step and because of the A-rule Req.® and of the A-rule effect Rcm.(|3]) the new extreme 
will be of magnitude bounded by 

El = Aai 

where full superscript ^ is used since the relocation has taken place. 
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Figure 4: The resulting situation at the head of the shock at the end of the 1-st time step. Only 
one extreme exists in this time step and it is of magnitude El . The numerical solution -before the 
remcshing procedure takes place- is depicted red, the new nodes -that occur after the remeshing 
procedure- are depicted in blue. The new node -that is closest to the old extreme- avoids the 
extreme by the A-rule resulting in a magnitude bounded by El < Xai. The reconstruction 
of the numerical solution over the new mesh results in the clipping of the magnitude of the 
extreme according to the A-rule. 



Remark 6. Regarding the continuation; the 1-st extreme shall "pollute" its neighbour by 
provoking the appearance of a 2-nd extreme of the opposite direction. 

Remark 7. (E^ Notation) We denote by E^^^^^ (half superscript) the bound on the magnitude 
of the m-th extreme at the k-th time step after time evolution and before the mesh reconstruction 
procedure. We denote by E^^^ (full superscript) the bound on the magnitude of m-th extreme 
at the end of the k-th time step, that is after the time evolution and the mesh reconstruction 
procedure. 

2-nd step We refer Figure [3] for the situation at the end of the 1-st step. At the end of the previous 
step, we had one extreme of magnitude bounded by El = Xai. Due to the time evolution 
we expect the 1-st extreme to evolve to a new value. We also expect the creation of a 
2-nd extreme at the left side of the 1-st extreme. We will study each extreme separately. 

1-st Extr. The Evolution Req.((T]) dictates that this extreme shall evolve as: 

Wl"-'^' - ul\ < Cmaxll^i - {^,,1 \ul ~ ul,\} 

where we use half superscript in u since the relocation procedure has not taken place 
yet. From the previous time step we have that, 

\ul - ul-i\ < El and - ul+, \ < 2EI + 02 

To justify the second inequality we return at the end of the time step k ~ 1 and 
notice that the node i -I- 1 is placed along the shock, which is -by symmetry- of 
variation at most El -\- TV{u°) -\- El. So the Evolution Req.((T]) for the 1-st Extreme 
reads, 

1^1+1/2 _ yij < C{2El + 02) = 2CE\ + 02 
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where we have defined 02 = Ca2. If now we set v to be the level from which we 
measure the magnitudes of the extremes (in this case the top of the shock), the 
previous bound recasts, 

\{u]+'^'-v)~{ul~v)\<2CEl + a2 

By setting E^^^^'^ = u^^^^^ — v and since El = uj —v we deduce that the magnitude 
of the 1-st extreme will be bounded as 



^1+1/2 ^ ^1 ^ 2C£;i 



02 



Now the relocation procedure takes place and because of the A-rule Req. ^ and of 
the A-rule effect Rem. ^ , the magnitude of the 1-st extreme at the end of this step 
shall be bounded as follows, 

E^ ^ X{El + 2CEl+a2) 
2-nd Extr. The Evolution Req.([T]) dictates that this extreme shall evolve as, 

- w^-il < C max - u-^sL " 1} 

where again half superscript is used on u^^^^'^ since the relocation procedure has not 
taken place yet. From the previous time step we know that 

\ul^i ~ ul-2\ = and \ul_^ - ul\ < El 
So the Evolution Rcq.([T]) recasts, for the 2-nd extreme as follows, 

\ul^l^'-uU<CEl 

or by noting that ul_i = f is the level from which we measure the magnitudes of 
the extremes, the previous bound recasts 

eI+'/' < CEl 

where, as we explained earlier half superscript is used because the relocation proce- 
dure has not taken place yet. 

Now the relocation procedure takes place and the A-rule Req. ([21) dictates that the 
magnitude of the 2-nd extreme at the end of this step shall be bounded as follows, 

Ei = XCEl 

So at the end of the 2-nd step the bounds on the existing extremes are as follows, 

Ef = X{El + 2CEl + aa), E'i ^ XCEl 
The situation at the end of this step is depicted in Figure [H 

Remark 8. The 2-nd extreme shall provoke the appearance of a new extreme of the 
opposite direction. This is the pollution process. 

3-rd step At the end of the previous step we had two extremes with magnitudes bounded by Ef 
and E2, sec Figure [S] In this step we expect them to evolve to new values Ef and we 
also expect a new extreme to appear, namely 
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Figure 5: The resulting situation at the end of the 2-nd time step. Two extremes of magnitudes 
E1 and El exist in this time step. The numerical solution before the remeshing procedure is 
depicted in red and the new nodes -after the remeshing procedure are depicted in blue. The 
new nodes -that are closest to the previous extremes- avoid the extremes by the A-rule 



1-st Extr. Following the discussion of the previous steps, we note that the evolution of the 1-st 
extreme will be governed by the Evolution Req.(IT]), so 

\u]^"' ^ul\<Cn.^^{\ul-uUm^uU\] 
where from the previous time steps we note 

\n^-uti\< El + El and \ul - ul^, \ < 2EI + 

We also note that from the previous time step the bound El = ^{El + 2CE\ -\- 02) 
is obviously larger than the bound El = XCEl hence the Evolution Req.([T]) for the 
1-st extreme reads as follows (after the subtraction of v), 

^2+1/2 ^ ^2 _^ 2CEl + as 

Now the mesh reconstruction procedure takes place and the new magnitude of the 
1-st extreme shall be bounded by 

Ef ^ X{El + 2CEl + as) 

2-nd Extr. Using similar arguments as before, the Evolution Req.((T]) dictates the evolution of 
the 2-nd extreme as follows, 

I 2+1/2 ^2 I ^ ri-2 -2 I 1-2 -2n 
K^i - w^-il < max||w,_i -Ui_2|, - ujj 

From the previous time step we note 

- "?-2 1 < El and \ul_, -ul\<El + El 

hence the Evolution Req.lJT]) for the 2-nd extreme yields, 

E^-^'/'- < El + C{El + El) 
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Figure 6: The resulting situation at the end of the 3-rd time step. Three extremes of magnitude 
El, i?2, E^ exist in this time step. 



Now, relocation takes place and the new 2-nd extreme shall be of magnitude bounded 

by 

El = XiEl + C{Ei+Ef)) 

3-rd Extr. Repeating the work we did for the 2-nd extreme in the previous time step, the 
magnitude of the 3-rd extreme after both the time evolution and the node relocation 
procedure will be bounded 

E^ = XCE2 

So at the end of the 3-rd step the bounds on the magnitudes of the existing extremes are 
as follows, 

Ef ^ X{Ef + 2CEf + as) = X\l + 2Cfai + A2(1 + 2C)a2 + Aag, 
El = X{El + C{El + El)) = X^2C{1 + 2C)ai + A^Caj, 
E^ — XCE2 — X'^C'^cLi 

Figure El depicts the situation at the end of the 3-rd step. 

For the sake of completeness we define the increases a.i that we used throughout the previous 
paragraph. For this we first analyse the variation of the shock at the fc-th time step. It consists 
of three parts, the oscillatory part at the top of the shock with magnitude at most i?f , the 
main part of the shock which is of variation TV{u'^) and the oscillatory part at the foot of the 
shock being of magnitude at most E^. 

Definition 3.1 (Definition of the increases). Let u'l he the value at the top of the shock. 
The node x^j^^ is located along the shock, in one of the three parts that consist the shock. 
We define 

«fc-(l«f-"ml-2i?i') + 
where the subscript + denotes the positive part. Moreover we define = Ca^. 

Remark 9. By definition, Ok describes the possibly more that 2£'^' distance — ''^i+il- That 
is, if — ""i+il < 2£'f then ~ 0; hence — 0. 
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We can generalise the situation in the k-th time step as foUows, 

/c-th step To generahse the description that we presented in the previous time steps, we introduce 
here the recursive relations for the general extreme m at the general time step A; > 1, 

= A (S^i + C ■ {E^-' + E^-_\)) , for m > 1 
E'i ^ X {E^-^ + 2C ■ E^-^ + flfc) , for m = 1 ^ ^ 

Remark 10. From the second equation if Rel. 0), it is evident that we add at least 
2CEI^^ in the increase of the 1-st extreme even if the actual increase needed for the 
highest node is less. This increases the magnitude of the 1-st extreme hut at the same 
time simplifies the presentation and the route of the proof. More precise increases, result 
in sharper final hounds. 

In analysing these recursive relations, we see that for the evolution of the extreme m 
-from the value E^^ to the value we take into account the neighbouring extreme in 
the right hand side, E^~\. To justify such a choice, we have to prove that the bounds 
E^ on the magnitudes of the extremes constitute a decreasing sequence with respect to 
TO = 1, . . . for every step k. That is we need to prove that Ei > > E^ > ■ ■ ■ . This is 
accomplished by the following lemma, 

Lemma 3.1. For every step k the magnitudes of the hounds given hy Rel-lBj) are in a 
decreasing order w.r.t m ~ 1, . . . 

Proof. Let's assume that in the step k, < E'^-^ for every m = 1, • • • . We shall prove 

that < E^^ for every to = 1, . . . 

TO = 1 The recursive relations Rel.lO read for E^^^ and £^2^^ as follows, 

E^+^ = \{E^ + CS^ + CE'i), 
eI+^ = X{El + CE'i + CE^ + flfe+i) 

Utihsing the induction hypotheses and that au+i > the result £'2^^ < E^^^ is 
immediate. 

TO > 1 The recursive relations Rel.® read for E^:^^ and E^^^ as follows, 

Et+'^\{Et+CEi^^ + CEt_,) 

The induction hypotheses states that < E^n — -^m.-ii ^'^ immediately we 

conclude that £;,^+_\ < E^^. 

It is so proven that in order to bound the new magnitude of every extreme we could use 
the recursive relations Rel.(|n]). □ 

Having devised recursive relations for the evolution of the extremes, that is Rel.® we 
continue with the study of the bounds of their magnitudes. 
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3.2 Extremes 

In this paragraph wc solve the recursive relation ReLlO, for every extreme m and we provide 
uniform -with respect to the time step k- bounds on the magnitude of the extremes. 

We start by providing bounds on the extremes of E^^ with respect to the sequence of increases 

Lemma 3.2 (Magnitude of the 1-st extreme). The magnitude of the first extreme in the k-th 
time step is bounded by, 

k 



a J, 



or, by setting I = k — j 



E^ < A ^ A'(l + 2C)'afc_i (6) 

Proof. By induction. We note from the previous discussion that 

1 

El < Xai = A^ Ai-J (l + 2C)^-^aj 

For the induction hypothesis we assume that the magnitude of the 1-st extreme is bounded in 
the fc-th time step as 

fc-i 

E'l < A^A'(l + 2C)'afc_i 

Using the evolution relation ([5]) of the 1-st extreme, that is 

Sf+i <\{E'l + 2CE'l+ak+i) 
we can bound its magnitude in the fc + 1 time step, 

E'-+^ < A ((1 + 2C)E'[ + ak+i) 
The right hand side recast -by the induction hypothesis- as follows, 

E'l+^ < A I (1 + 2C)A^ A'^-^ a + 2C)^-^aj + ak+i 

< A A^-+i-^(l + 2C)'^-+i--'a, + A^+i-f'^+i'll + 2C)''+^~^''+^'>ak+i 

k+l 

< A^A'^-+i-^'(l + 2C)'^+i- 

This completes the proof regarding the bound of the magnitude of the 1-st extreme. □ 
Wc now need a similar bound on the magnitude of the 2-nd extreme, 
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Lemma 3.3 (Magnitude of the 2-nd extreme). The magnitude of the second extreme in the 
k-th time step is bounded by, 



or by setting I = k — j 

< \'CJ2 (i [ ^\'-\l + 2Cy-'au-l (7) 
1=1 ^ ^ 

Proof. Proof by induction using relations dS]), dH) and the fact (^) + (^."J = (";J:J). We 
note from the previous discussion that 



For the induction hypothesis we assume that 

fc-i 



i=i 

and for the induction step we have the following 

E^^'^^ = A [El^ + C(E^ + E'l)) = A ((1 + C)E^ + CE'[) 

Where in the last step we utilised the induction hypothesis. Now, since 1 + C < 1 + 2C the 
bound recasts 

k~l / ^ X k-1 



E^+' <xixcJ2 jA'(l + 2C)'afc_z + AC^A'(l + 2C)'afc_z 



(E ((; ! 1) + Q) ^'(1 + 2C)'«.-. + a. j 
A^C ("xi ) ^) A'(l + 2C)'a,_, + (° -^"^^ + 2C)"a.._o' 



fe-i 



/ + 1 



A^C^C ; ^)a'(1 + 2C)W-; 
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Finally we set fj. = I + 1 and the bound on the magnitude of the 2-nd extreme reads, 

(fe+i)-i 



(8) 



and this completes the proof regarding the magnitude of the 2-nd extreme. □ 
Similarly we prove that the magnitude of the 3-rd extreme k-th time step is bounded by, 



or by setting I = k — j, 

fc-i 

J -2 



1=2 



We can generalise the previous lemmas, in a compact form for the m-th extreme in the fc-th 
time step. 

Lemma 3.4 (Magnitude of the m-th extreme). The magnitude of the m-th extreme in the k-th 
time step is hounded by, 

Et < A™C""-i ^ (fc - j-L + l) ^'^'"^^"^^1 + 2C)'^'~^-™+^aj, 
or by setting I = k — j , 

Et < A"C"'-i ^ (l- !n l) ^'""^'(1 + 2C)'-"+iafc^i (9) 

l=m-l ^ ™ ^ ^ 

Proof. The proof is exactly the same as in the 2-nd extreme and is omitted. □ 

The last lemmas provided bounds on the magnitudes of the extremes. In the following 
remark we merge these bounds in a single relation valid for every m = 1,2,3,.... 

Remark 11. The bound we have extracted for the m-th extreme at the k-th time step, that is 
Rel.^: 

£;^<A'"C""-i j^^_^')A'-™+i(l + 2C)' 



l—m—1 



is valid for every m ~ 1, 2, 3, . . -not just for m > 1- since, for m = 1 the bound we extracted 
for the 1-st extreme at the k-th time step, that is Rel.(^, 



fe-i 



^i' < A^A'(l + 2C)'afc_, 



1=0 

can be written in the form 

fc-i 



1=1-1 



E'l<X'C'-' Y (^_|^ jA'-i+i(l + 2C)'-i+ia.._, 
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So far we have described the creation and evolution of the extremes. We provided Recursive 
relations ([5]) that connect the magnitudes of extremes, we have solved the recursions and merged 
the magnitudes of the extremes into a single relation Q. 

Remark 12. The bound we have extracted for the m-th extreme at the k-th time step, that is 
Rel.(^: 

m-l , . s 



is valid for every m ~ 1, . . . since the bound we extracted for the 1-st extreme at the k-th time 
step, that is Rel. ^B^: 



since the former can be written in the form of relation (0) for m = 1 . 

Now, we note that the bounds on the magnitudes of the extremes, that we have extracted, 
depend on the time step k. We shall bound the magnitudes of the extremes uniformly with 
respect to the time step k. This will allow us to provide the final proof regarding the total 
variation increase due to oscillations. 

Lemma 3.5 (Uniform -with respect to the time step k- bound on the extremes). // there is 
a constant M > such that ai < CM for every i = 0, ■ ■ ■ oo and if X -\- 2XC < 1 then every 
extreme m is uniformly -with respect to the time step k- bounded as, 

Proof. The magnitude of the bound of the m-th extreme, m = 1, 2, 3, ... at the fc-tli time step, 
with m < fc, is given by the Rel.©, 

Since the increase Oi are uniformly bounded, a,; < CM (we refer to the definition of the increases 
Oi Def. (|3.ip ) we can bound the extremes as, 



Et < y^C^M V , A'-™+Vl + 2C) 

^— ^ \ / — 771-1-1 / 



A — 771+1 

l—m—1 

Setting u = I — m + 1 the previous relation reads 

k- 



\/ — m+1 



Et < A"C"M ^ J ^ + ] A'^(l + 2Cy 



u=0 



All the terms inside the sum arc positive, hence the right hand side of the previous relation is 
increasing with respect to fc. Hence it can be bounded for fc = cxo as follows. 



v -\- m — 1 



i/=0 



£;^<A™C"M^ A''(1 + 2C)'' 
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or ^ 

For the convergence of the previous infinite sum we recall at this point the power series expan- 
sion. 

Efv + m — 1\ ,, 1 , II-, 

]t = ^ — , whenever \t\ < 1, 

and since A + 2AC < 1 the last bound on E!^ reads as follows, 



= M 



(1 - A - 2\cy 

AC 



1 - A - 2\C , 

Which proves the assertion of the lemma. □ 

Remark 13. // moreover we assume A + SAC < 1 -instead of X + 2AC < 1- then the sequence 
of bounds on the extremes {E^} is decreasing with respect to m, since they can be written as, 



E^, < M 

hence 



AC 



1 - A - 2AC 



lim Et < lim M f - ] = 



1 - A - 2AC, 

•xe E^^ > and the fraction i_x-2\c ^ ^ because A + SAC < 1. 
We are ready now to measure the total variation increase due to the oscillations. 



3.3 Variation 

In the previous lemma we proved that each extreme separately is of bounded magnitude, uni- 
formly with respect to the time steps k. The next theorem is the basic one and states that 
in addition to the magnitude of the extremes, also the sum of the extremes is also bounded 
uniformly with respect to the time step k. 

Theorem 3.1 (Main Result). We assume that the requirements Req.(^ and i?eg. (0j are sat- 
isfied for A such that A -I- SAC < 1 by the numerical scheme and the mesh. We more over 
assume that the sequence {ai,i = l,oo} is uniformly bounded Oi < CM. Then the sum of the 
magnitudes of the extremes is uniformly -with respect to the time step k- bounded as follows, 

E ^™ ^ ^ ^-A-SAC 

m— 1 



Proof. We shall utilise relation (|10p . which is valid since the requirements of the relevant lemma 



arc satisfied. At the end of the fc-th step we have k extremes E^, . . . , E^. The sum -with 



respect to m- of their magnitudes can be bounded as 

fe fc 



E 



rn—1 m— 1 ^ ^ rn—1 ^ ^ i — A — zAO 

, 1 - A-2AC 

< Af 

- 1 - A - SAC 
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where the second inequahty and the equahty are vahd since A + SAC < 1. □ 



Summarising and concluding we can state the following theorem, which constitute our target 
result on the Total Variation Increase. 

Theorem 3.2 (Total Variation Increase Bound). Given the requirements of the previous The- 
orem, the Total Variation increase due to the oscillations is bounded and given by 

Proof. The variation of the oscillatory part is bounded by twice the magnitude of the extremes. 
So, 

AC \" „. 1-A-2AC 



TVI < 2A/ ^ (^-^^ j <2M 



1 - A - SAC 

where the last inequality results from the previous Theorem (Main Result). □ 



Although we proved our main result, we can gain better insight if we study the contribution 
each increase factor a,; has on the total variation. For this reason we include the following 
paragraph. 



3.4 Variation- Revisited 

We shall follow now another approach that will provide us a with further insight of the " pollu- 
tion" process and with a sharper bound on the Total Variation Increase. 

This approach differs from the previous one in the sense that instead of adding directly 
the magnitudes of the extremes E^-^, we compute the contributions of the increase terms a;, 
z = 1, 2, 3, ... in the each one of the extremes E^-^ separately. Then we add this contributions 
with respect to ai. 



ai cont. The contribution of ai in the fc-th step. 

In the fc-th time step there exist k extremes and the increase factor ai is present in each 
one of these extremes. So we extract the contribution of ai in all the extremes that are 
produced during this procedure up to the fc-tli time step. 

The contribution of oi in the 1-st extreme in the fc-th time step is given by the relation 
([6]) and reads as 

and in the general extreme m the contribution of ai is 

\k — m J 

Summing these contributions with respect to to we end up with the total contribution of 
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ai in the k-th time step, 

jk 



(for V = k — m) 

cont. The contribution of a2 in the fc-th step, 

Similarly we notice that in the fc-th time step the increase factor 02 contributes in all the 
extremes except the last one, m ^ k and its total contributions is 

= A(A + 3AC)'=-2^2 

cont. The contribution of a„i in the k-th step. 

More generally, in the fc-th time step the increase factor a„i contributes in all but m — 1 
extremes (the last m — 1) in the fc-th time step (fc > to), and its total contribution is 

I^^ = A(A + 3AC)'=-"a„ 



V fc — TO / 

m=l ^ 

= A^ V f ^" " M C"-' (1 + 2C)'=-™ai 

m=l ~ 

= A''(l + 3C)''-iai = A(A + 3XC)''-^ai 



Remark 14. FKe note here that as long as A < i.e X + SAC < 1 each one of these 

contributions converges to as fc ^ 00, 

lim = lim A(A + 3AC)'=-'"a,„ = if A < — ^ (12) 

T/iis is the very essence of the X-rule effect. Namely a remeshing procedure which respects the 
X-rule requirement (0) at the extremes is able to limit the increase of the variation due to each 
Oi -eventually kill it- and hence provide us with a control over the total variation of the scheme. 

To finalise this second approach to the Total Variation Increase due to oscillations we con- 
tinue by summing the contributions of all the a^'s in the fc-th time step. This will result in half 
the Total Variation Increase due to oscillations in the fc-th time step. 

By the previous talk the following corollary is obvious. 

Corollary 3.1 (Total contribution in th fc-th step). In the k-th time step we have contribution 
by ai, 02, . . . , Ok, with sum, 

k fc 

IL = E = ^ E + 3AC)'=-™a™ (13) 

Tn=X m—1 

Corollary 3.2 (Result 1). If we assume that the sequence Oi, is bounded i.e there exists M > 
such that ai < CM for aZH = 1, . . . , 00 and that X + 2AC < 1 then the total contribution in the 
k-th time step reads, 
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Proof. Since the increase factors Ui are uniformly bounded as < CM, their total contribution 
given in Rel. p^ reads 



fc-i 

= fe-T 



< XCM E + 3AC)^-" " = " = ACA/ ^(A + ZXCf 

rn— 1 

1 - (A + 3AC)'^ 



n=0 

I I XI — t— . I I / ^ ^ 

XCM- 



1 - (A + SAC) 

the previous sequence, in the right hand side, is increasing with respect to the time step fc, so 
by taking the limit as fc — > cx) we deduce an uniform -with respect to k- bound on the total 
contribution 

= Hm II, < XCM hm ^ ' (-^ + ^AC)' ^ AC 

*°* k-,oo *°* - fc^oo 1 - (A + SAC) 1 - (A + SAC) 

This is exactly the result we were looking for since now, the Total Variation Increase due to 
oscillations is bounded 

2AC 

TVT < 9 ■ T°° < M 

iVi<^ i_(A + SAC)^^^ 

□ 

Corollary 3.3 (Result 2). If we assume that the sequence ai, i = 1 ... oo is uniformly bounded 
as ai < CM = CTV(uo) then the previous bound on the total variation increase becomes, 

Remark 15. This result is even better than the previous one given Rel. ^4^ since the bound 
that provides on the increase of the Total Variation due to oscillations is directly related to the 
variation of the initial condition uq. 

With more delicate assumptions on the increases ai we get the following Corollary, 

Corollary 3.4 (Results). If we assume that the sum of all increases Oi is finite '^'i^^Oi — 
A < oo then the Total Variation Increase due to oscillations diminishes with respect to the time 
step k. 

Proof. We note that the following sums converge 

oo _j 

i:(^ + ^^^)'- l-(A + SAC) <- 



a,; < oo 

i=0 



Moreover the sum 



oo oo I k \ oo I k 

E ^tot = E hE(^ + 3AC)^-^a, = Aj] K](A + SAC)'=--'a, 
k=i k=i \ j=i J k=i \j=i 
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constitutes the sum of the terms of the Cauchy product of the series X^i^i '^i ^^'^ ^iLii^ + 
3XCy. So we deduce that the foUowing sum also converges, 



k=l 



tot < OO 



and hence the sequence I^^f must converge to 0, so 



^ hm = 

A; — 'OO 



□ 



Remark 16. The last result states that if the sum ^ is finite then the overall increase of the 
variation due to the oscillations produced by the ais' diminishes with respect to the time step k. 

Moreover we note that the requirement 'Y^ai < oo is supported numerically, since the first 
node on the shock at the right hand side of the first, is always very close to the first extreme. 
This is due to the high density of the mesh at the shock. Hence by the definition of the increase 
factors Def. H3.1\) . Ok = and so Ok = 0. 

We have devised two different bounds concerning the Total Variation Increase due to oscilla- 
tions. The first was immediate summation of the bounds of the magnitudes of the extremes 
and resulted in the bound Rel. (fTT|) . The second one came by investigating the contributions of 
the increase factors and resulted in the bound Rel.flT 



Comparison of the two bounds We start this paragraph by restating the two bounds on 
the Total Variation. The first given in Rel. (fTT|) 



and the 2-nd given in Rel.([T 



1 - A - 3AC 



B. ^ M 
1 - A - SAC 



Before the comparison, a comment on the nature of the second bound 
To compare the two bounds one can examine their ratio, that is the fraction of the bound 
RcLdHI) over the bound of Rel.((TT]), 



Lemma 3.6. (Comparison of the bounds B2 and Bi) If X + 'SXC < 1 then ^ < 1. // 
A + 4AC < 1 then |f < 5 

Proof. The ratio of the bounds B2 and Bi is 



moreover 



„ 2ACM ,^ 

-D2 1-A-3AC '^L/ 



2A/i5Ml§ 1 - A - 2AC 

If A + SAC < 1 then ^ < 1 since AC < 1 — A — 2AC, hence the later bound B2 is sharper. If 
moreover A + 4AC < 1 instead of A + SAC < 1 then # < 1 since 2AC < 1 - A - 2AC. □ 

The meaning of this proposition is that by a careful selection of the respect factor A the 
bounds on the increase of the variation in the second approach can be significantly better than 
that of the first approach. 
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4 Main Adaptive Scheme (MAS) 



In this section the procedure that is responsible for the creation and manipulation of non- 
uniform meshes, their connection with the numerical solutions and the time evolution. We 
start by providing the definition of MAS once more, 

give references 

Definition 4.1 (Main Adaptive Scheme (MAS)). Given, at time step n, the mesh M" = {a = 
Xi<---< = b} and the approximations U" ~ {^i i ■ • ■ i ^n}' steps of the (MAS) are as 
follows: 

1. (Mesh Reconstruction) 

Construct new mesh M^+i = {a = < • ■ • < a;^+^ = 6} 

2. (Solution Update) 

Using the old mesh M" the approximations C/" and the new mesh Af"+^; 

2a. construct a piecewise linear function V"'{x) such that V''^\m^ = 

2b. define the updated approximations C/" = {m", . . . ,u%} as U" = ^"l^^+i 

3. (Time Evolution) 

Use the new mesh Af"+^, the new approximations t/" and the numerical scheme to march 
in time and compute U^'^^ ~ {u^^^, . . . , u^^} 

4. (Loop) 

Repeat the Step 1.-3. with MJ^'^^ , U"^^ as initial data. 

The Mesh reconstruction procedure (Step 1.) is a way of relocating the nodes of the mesh 
according to the geometric information contained in a discrete function. The basic idea of the 
mesh reconstruction is simple and geometric 

"in areas where the numerical solution is smoother /flatter we need less nodes where, in 
contrary, in areas where the numerical solution is less smooth/fiat more nodes are in order" 

The key point in this procedure is the way that we measure the geometric information of 
the discrete function. This is accomplished by using two auxiliary functions. The first one - 
estimator function- measures the geometric information of the discrete function and the second 
one - monitor function- redistributes the nodes according to the information measured by the 
estimator function. 

Some examples of estimator functions are the arclength estimator, the gradient estimator 
and the curvature estimator. Throughout this work we shall use the curvature estimator. 



The curvature estimator function To start with we consider a smooth (in the classical 
sense) function U. The function Ku{x) that measures the curvature of the smooth function U 
is defined as follows, 

\U"{x)\ 



Ku{x) 



{l + {U'{x)Yf^ 



The functions that we deal with (numerical solutions) are not smooth, they are merely point 
values and hence a discrete analog of the curvature estimator functions is needed. To gain a 
discrete analog of the estimator function we just have to discretizc the derivatives that appear 
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in the smooth estimator. Again several choices are possible, but the one that we use throughout 
this work is the following, 




where MJ^ — {x", i — 1, . . . , N} and [/" = {u", i ~ 1, . . . , iV} are the non- uniform mesh, and 
the approximate solution at time step n. 

With the discrete estimator we can evaluate the respective geometric property of the numer- 
ical solution itself. By performing a point by point evaluation of this discrete estimator results 
in a finite sequence that contains the measured information of every node (x", u"). That is. 

Remark 17. We note here that the proposed mesh reconstruction procedure, assumes positive 
values on behalf of the discrete estimator. So now on we assume that Kf^"^^ > -even though 
this is obvious for our discrete curvature estimator. 

Remark 18. In areas where the numerical solution is flat, all the nodes yield Kf"'^'^ = 0, 
hence no information is extracted. To avoid such a case we select a e > and set Kf'^'^^ = 
m.ax{Kf^''^,e} for every i. A typical value of e k, 10~^^. 

Remark 19. In the initial condition of a Riemann problem all the nodes ( except for the 2 nodes 
at the top and bottom of the discontinuity) attain the minimum "information" Kf'^'^'^ = e. The 
other two nodes attain very large "information" . This abrupt change of "information" yields 
a non smooth non-uniform mesh. To avoid such abrupt information change we use a constant 
pw > and set Kf'^'^'^ = (^xf^"')^^ for every i. This results in a smoother transition of 
consequent discrete estimator values, a typically value of pw ~ 0.9. 

Finally, these points, i.e {(x", Kf'"^'^), • • • , (x^, A'^'*'^'')} are interpolated by a piecewise lin- 
ear function I^dsar. 

We move on to the monitor function that will provide us with the new nodes. 



The monitor function The evaluation of the monitor function starts from its discrete analog, 
that is we first evaluate the monitor function in every old node xf and the we construct the 
continuous monitor function by linear interpolation of the discrete monitor values. 

We integrate the piecewise linear function Ij^dscr to find the value of the discrete monitor 
function in every node x", 

= I ' Iu^[x)dx. 
Jo 

This results in a new sequence, 

{ix^,M:l)a = 0,--- ,N)] 

This sequence is positive and strictly increasing since Kf'^'^^ > for every z = 0, ■ • ■ , N. 

Finally, we interpolate over the values of this sequence by a piecewise linear function Mijn (x), 
which is continuous, positive and strictly increasing and so attains it maximum at the right end 
Af(l). We are ready now to relocate the nodes of the mesh. 



22 



0.6 0.8 



0.48 0.485 0.49 0495 0.5 0.505 0.51 0.515 0.52 



Figure 7: A typical estimator function is depicted along with the respective numerical solution. 
The right graph is a focused version of the left one. 
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Figure 8: A typical monitor function is depicted along with the respective numerical solution. 
The right graph is a focused version of the left one. 
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Figure 9: We see in this graph the result of the relocation procedure. The initial function 
before the mesh reconstruction procedure takes place is depicted in red. In blue we can see the 
reconstructed numerical solution, after the mesh reconstruction procedure. We can see that the 
density of the nodes is higher around the area of interest of the numerical solution. We also 
notice that there are nodes placed along the slope of the shock. 
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Mesh reconstruction What wc now need is a new set of nodes, {x"^^, i = 0, • • • , N}, with 
Xq^^ = such that they equi-distribute the total information M(l) that we measured. This is 
accomphshed by solving -recursively with respect to a;"_J*]^- the system, 

x^+^ =0, 

M,,„(.<+i)-A%„(xri) = iA.%„(l) ' 

Remark 20. It is obvious the last new node x^^ shall be the right end of our interval, that is 
x^^ ~ 1. Moreover the fact that Mijn{x) is strictly increasing, allows for inversion and since 
it is piecewise linear reduces the computational cost of the solution of the previous system to 
0{N). 

Remark 21. Let us note that the number N of nodes is constant. Let us also note that the 
mesh reconstruction procedure is not related with the numerical scheme that we use for the 
evolution part of the problem, it is merely connected to the geometry of the numerical solution 
itself. 

In Figures [7] and [S] a typical Estimator and Monitor function are depicted respectively. In 
Figure [H one can see the affect that the node relocation procedure has on initial condition. 



5 Computational considerations 

In this section we describe the numerical implementation of the Coupling Req.(l3]), and the 
connection of the A-rule Reg.© with the A-rule effect -as described in Rem.Q- when several 
nodes are located in between two consequent extremes. 



Numerical implementation of the coupling requirement We start by restating the 
coupling requirement, 

Requirement (Coupling requirement). The constants C of the evolution Req.^^ and A of the 

mesh reconstruction Req. ^ are connected via the following relation: 

A + SAC < 1 (15) 

In this section we shall discuss the numerical implementation of this requirement. We start 
by defining /" to be the set of indices of the nodes a;"^^ that are placed -after the mesh 
reconstruction step- "close" to a position of a local extreme of f7", that is 

/" = {j I x^^^ e [a;",x"+i) for some i and J7" exhibits local extreme at cc" or 

For every j £ ?7" exhibits a local extreme either at xf or xf^^, so we set: 

^. = ^^(1 + 3C) 



or respectively 

A^-'^n ^(1 + 3C) 



where the constant C is related to the numerical scheme under discussion. With this notation 
the coupling requirement reads for the discrete case as follows: 

max A, < 1 
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Figure 10: A typical evolution of the maximum and the average Aj, j £ /" with respect to the 
time steps n (horizontal axis). In the left graph we present the full range of their evolution, 
while the right graph is a focused version of the left one where one can more clearly see that 
the value 1 is an upper bound for the maxj-g/n Aj, the same is true for the average Aj, j & Ij'- 



To impose numerically this requirement we check, for every j G whether Aj > 1, and if 
so, we correct accordingly the position of the node x'^'^^- More specifically, if the local extreme 
is at the node a;" we set 

where < e < 1, with a typical value of e = 0.2. 

Similarly we treat the case where a local extreme is at the node xfj^i- We refer to Figure 
[TO] for a typical graph of the maximum and the average Aj, j S /" with respect to time step n. 



Connection of the A-rule requirement with the A-rule effect We now explain the 
connection of the A-rulc Rcq.(l2]) with the A-rulc effect -as described in Rem.(|4])- when several 
nodes are located in between two consequent extremes. 

For this, it is sufficient to show that an application of the mesh reconstruction (Step 1.) 
and the solution update (Step 2.) has a diffusive effect on the discrete numerical solution; since 
then, continuous repetitions of Step 1. and Step 2. -which do not consume any time of the 
physical evolution of the problem- yield the required A-rule effect Rem.dH). 

We refer to Figure [TT] and assume that a piecewise linear function attains a local maximum 
at the node with value u". The neighbouring nodes are x"_i and xf^i with values u^_i, 
■""-(-1 respectively. The slopes, left and right of the node x" arc given by 



and A^ 



where = x^ — xf_i and hf = xf_^-^ — xf. 

Next we assume that there is a smooth function u that interpolates the points {xf_i,uf_i), 
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Figure 11: The numerical solution before the mesh reconstruction is (x", w") and is depicted in 
red, after the mesh reconstruction is (x""''^, u" ) and is depicted in blue. The solution is updated 
by projection of the new nodes on piecewise linears. The projected value of the old node a;" on 
the linear new (blue) linear segment is u". 



(a;",M") and {x"^j^,uf^i), and we compute 

— _ hf hf _ hf hf 

h- + ht ~ 



K +ht K +h+ 

K + 

'-u"ix^)^9mi)±2wm 

2 ^ K+hf 



(16) 



We perform a mesh reconstruction and a solution update step which yield the new mesh 
M"+^ and the update solution t/". We define 

h^^x-~x-+\ h^^x^ll-x^ 

and note that due to the linearity and the interpolation the following is valid 



hi 



n 



We interpolate linearly over the points (x"^ , u"), {x^Jl^i , ^^j+i) smd compute the new value 
of the previous node (cf Figure (fTT|) ). To do so. we note that G {x"'^^ , x'^^^) and is given 
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by the following convex combination 



and so the projected value m" of the node over the new linear segment is given by the same 
convex combination with respect to and Wj+i, 



hi+h2 ' ^ hi+h2 ^+' 

" ■ -(Ar-Ar_i) (17) 



= K + ^u"{x2) + -^^^ mKf) + 0{{htf)) (using RcLdUl)) (18) 
where C = ^(V + ^+)>0. 

Remark 22. i?e/. {T^j reveals the diffusive property of the mesh reconstruction (Step 1.) and 
the solution update (Step 2.). We moreover note that in the case of a uniform mesh either 
hi = or h2 = 0, which both yield C — 0. 

As explained earlier, this remark is useful only if there are several nodes in between two 
consequent extremes. In this case, continuous repetitions of the mesh reconstruction and the 
solution update guarantee that the extremes will satisfy the A-rule effect. We note though that 
in practise, applying the mesh reconstruction procedure and the solution update once in every 
time step is sufficient. 



6 Numerical tests 

As stated at the introduction the numerical schemes that we shall discuss are oscillatory, either 
due to their dispersive or to their anti-diffusive nature. Here we shall only state their description 
for non-uniform meshes and prove that they satisfy the Evolution Requirement Req.(IT|), which 
we restate 

Requirement (Evolution requirement). There exists a constant C > independent of the time 
step n and the node i such that, 

- <| < Cmax {K+i - <|, |< - (19) 

The problems that we shall deal with, are the Transport Equation 

Mt -hux = 0, a; e [0,1], t e [0,1] 

and the inviscid Burgers Equation 




= 0, xe [0,1], t e [0,1] 



both with jump initial conditions 

uo{x) = A'[o,i/2](a;), x G [0, 1] 
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6.1 Richtmyer 2-step Lax-Wendroff 

In this approach we consider the non-uniform cell centered discretization of the domain in cells 



t^. = «-l/2,<+l/2) With \C2\=K 
The mesh MJ^ = G Z} consists of the middle points, 



X, = ■ — — hence x, - x, , = 



For this description of the grid we propose the following numerical scheme as the generalisation 
on non-uniform meshes of the Richtmyer 2-step Lax-Wendorff numerical scheme, 











i+l/2 













n+l 



At 



i^(/K-hi/2)-/("r-i/2)) 



Tt(/«+i)-/«)) 



or 



with 



i^.+ l/2 = /(^^^:Vl/2) = /( ,„+f5\„+l < + 



At 



Remark 23. A straight forward computation reveals that this scheme reduces to the usual 
Richtmyer 2-step Lax-Wendroff scheme when then mesh is uniform i.e /i" — h for every i,n. 

We need to bound the difference 



At 



|<+^-<|<^ /KVi/2)-/K-i/2) 



<CFL\u*^,/,-uU/,\ 



-1/2I 



and now the difference 



Fi+1/2 ~ ""1-1/2! 



■n+l 



At 



,ri+l , .n + l "i+1 7 n+l , .n+l (/("i+l) / ("i )) 



H+1 



' n+l I I,n+1 I 
h r "i+1 

h^+^ ^. 

htl + h'l+^ K-i + + 



I "i+i 
"-i-i 



At 



-(/«)-/«_!)) 



by setting = f-i and = M2 the previous bound recasts into: 
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Figure 12: In this graph are presented two time instances of the numerical solutions of the 
inviscid Burgers equation using the Richtmyer 2-step Lax Wendroff over uniform and non- 
uniform meshes. In the first row of graphs the full domain is presented, while in the second 
row the graphs are focused on the top of the shock. The uniform mesh case -depicted in 
blue- exhibits oscillations due to the dispersive nature of the numerical scheme where as the 
non-uniform mesh case -depicted in green- is clean. 



Ml< + (1-Ml)<+1 



- ^l2u7-l - (1 - ^^2K + 

<MiK-<+il+M2K-< 



- ifi^l+l) - /«)) 



J+1 
At 



"-i-i 



- (fiu-) - f{uU)) 



At 



H+l 



-max|/'||u^Yi 



At 



2 mm /i j 2 mm ti^ 

< (l + /ii +/i2 + CFL)max{|u"_i 

< (3 + CFL) max{|«r„i ~ vJ}\, - u^^,\} 



— max l/'llw" 



where the last inequality is valid since < /ii,/i2 < 1. So the overall boimd reads, 
- u7\ < CFL{3 + CFL) max { l^i^+i ~ ^i^l, K - uti\} 

The constant C in this case is chosen to be C = CFL{3 + CFL), for this choice the evolution 
requirement is satisfied. We refer to Figures [T^] and [T3] for comparative graphs between the 
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Figure 13: Transport equation with velocity a ~ I, using the Richtmyer 2-step Lax WendrofF. 
Again, oscillations appear in the uniform mesh case, whereas the non-uniform one is clean. 



uniform and non-uniform mesh case for the Richtmyer 2-step Lax- WendrofF scheme. 
6.2 MacCormack 

In this approach we consider the non-uniform cell centered discretization of the domain in cells 

cr = (xr_i/2,xr+i/2) with \ct\ = k 

The mesh M" = {xf,i G Z} consists of the middle points. 



For this description of the grid we propose the following scheme as the generalisation of the 
MacCormack, 
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Figure 14: Transport equation with velocity a ~ 1, using the MacCormack scheme over uniform 
and non-uniform meshes. Again oscillations are apparent in the uniform mesh case -due to the 
dispersive nature of the scheme- whereas the non-uniform is clean. 



We first rewrite the scheme in the following form, for /* = f{u*) and fi = f{uf) 



2 2 

~ T ^ ^ ^ 

~ ^ , n+1 , ;,n+l ^ ■^^) ^ I n+1 , (■^j ^ /i-l) 

So, to prove the Evolution Requirement for this scheme we need to bound 

Wi^ ^ I " I ^ ,n+l I I n+1 ^ fi) ~ .n+1 , .n+1 (/i ~ /j-l)l 

which reads, 

wr' -K\<^ (K+i - + 1< - <-ii) 

< (l^ir+i - K\ + \K - + CFL\u^ - u^ui + cfl\vj:+, 

< CFL(1 + CFL)max{|^ir+i - Kl \K " 
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Figure 15: Inviscid Burgers equation, using the unstable FTCS. The oscillations in the uniform 
case are due to the anti-diffusive nature of the numerical scheme. The non-uniform mesh case 
is clean. 

So, the constant C in this case is chosen to be C = CFL{1 + CFL) and for this choice the 
Evolution Requirement is satisfied. We refer to Figure FHl for a comparison graph between the 
uniform and non-uniform mesh case for the MacCormack scheme. 

6.3 Unstable Centered - FTCS 

This scheme produces oscillations due to its anti-diffusive nature, which property is also re- 
sponsible for the instability of the scheme. 

In this approach we consider the non-uniform mesh 

A4" = {x^ i e Z} with ^x'i- 
The middle points 2:"_j^2 — 2~'~^' define a partition of the domain in cells, 

c^r = «-i/2,<+i/2) with icri = ' ^ 

For this description of the grid we discuss the known to be unstable Forward in Time 
Centered in Space (FTCS) scheme 
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This scheme can be written in conservative form as follows, 



u 



n+l 



2At 



'H + "-i+l 



+ l(-^i+l/2 -^1-1/2) 



with 



^1+1/2 



Easily we deduce that 



K+'-u-\ < --— ^max|/'||<+i < —- (l^ir+i -<| + |< 

2 mm /ij ^ 

< CFL max { j^i^+i - "7 1 , 1"? ^ "^-i I } 

The constant C in this case is chosen to be C = CFL, for this choice the Evolution requirement 
is satisfied. We refer to Figure fTSl for a comparison graph between the uniform and non-uniform 
mesh case for the FTCS scheme. 



7 Conclusions 

In this work we investigated the creation and evolution of oscillations over non-uniform, adap- 
tively redefined meshes. The mesh reconstruction is driven by the geometry of the numerical 
solution and the solution update is performed by interpolation over picccwise linear functions. 
The numerical schemes considered are 3-point non-uniform versions of oscillatory numerical 
schemes. The overall process is dictated by the Main Adaptive Scheme (MAS). 

We prove under specific assumptions/requirements on the mesh reconstruction, the solution 
update and the numerical schemes that the numerical solution is of Bounded Total Variation; 
furthermore, under more strict assumptions, we prove that the increase of the Total Variation 
decreases with time. We also describe the diffusive property of the mesh reconstruction and 
solution update steps of MAS and finally we provide numerical test supporting the outcomes 
of this study. 
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